The Electronic Correlation Strength of Pu 
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A new electronic quantity, the correlation strength, is defined as a necessary step for understanding the prop- 
erties and trends in strongly correlated electronic materials. As a test case, this is applied to the different phases 
of elemental Pu. Within the GW approximation we have surprisingly found a "universal" scaling relationship, 
where the /-electron bandwidth reduction due to correlation effects is shown to depend only upon the local 
density approximation (LDA) bandwidth and is otherwise independent of crystal structure and lattice constant. 

PACS numbers: 71.10.-w, 71.27.+a, 71.30.Mb 



Many technologically important materials have strong 
electron-electron correlation effects. They exhibit large 
anomalies in their physical properties when compared with 
materials that do not have these effects, and have significant 
deviations in their electronic-structure from that predicted by 
conventional band-structure theory based on the local-density 
approximation (LDA). Because these anomalies and devia- 
tions are caused by electronic correlation effects, which often 
dominate the physics of these materials, in this paper we de- 
fine a new quantity, which we call "correlation strength," or 
C, as a necessary step in order to be able to describe trends 
and bring order into our understanding of these materials. We 
emphasize the word "quantity" since a quantitative measure 
is needed to answer the question: "how strong are the elec- 
tronic correlations?" Without some understanding of how big 
this is, it is not possible to make sense of the properties of 
these materials. In this context, "correlation" is defined in a 
way somewhat different from how is sometimes used (e.g., in 
the term "exchange-correlation potential"). By "correlation" 
we specifically mean "correlation beyond LDA theory". This 
usage reflects the way the term is often loosely used in com- 
mon terminology in the area of strongly correlated electronic 
systems. 

To create a new quantity requires determining a "scale" by 
which to measure its size. In principle, any experimental or 
theoretical property (e.g., specific heat) that monotonically 
increases or decreases over the full range of correlation ef- 
fects, where we define correlation strength to lie between 
for none and 1 for full correlation, can be used as a measure 
of this quantity. Hence correlation strength is an indetermi- 
nant quantity and depends on the property used to define it. 
However, this does not matter since only relative rather than 
any absolute strength is important for characterizing these ma- 
terials and for predicting trends in their properties. Any mea- 
sure based on one property can easily be converted to that 
based on another property. In this paper we develop a the- 
oretical correlation strength based on the GW approximation 
1 1 , 2 ] to electronic-structure theory and apply it to plutonium 
|3l|4[], which is known to have significant correlation effects. 
The GW approximation is named for the correction term in 



this theory, which is a Green's function G times a screened 
Coulomb interaction W. We also demonstrate a scaling rela- 
tionship that is universal in that it is independent of crystal 
structure and atomic volume. The ideas in this paper could 
certainly be modified and generalized to be able to treat other 
types of correlated materials (e.g., spin-fluctuation or high- 
temperature superconducting materials) by using other elec- 
tronic properties to determine a correlation strength and by 
using more sophisticated theoretical techniques than are con- 
sidered here. 

Our meaning of correlation makes it necessary to use a 
theory that includes correlation effects that go beyond those 
included by the LDA approximation in order to determine a 
theoretical correlation strength. This is challenging, since the 
most sophisticated treatments of correlation effects are mainly 
confined to abstract theoretical models. These parameterize 
the electronic structure in such an oversimplified manner that 
the connection with actual materials examined experimentally 
can be somewhat vague. Even recent methods, such as dy- 
namical mean-field theory (DMFT) still require fitting 
Hubbard parameters fl, and the model part, which is ad hoc 
and cannot be precisely defined or derived, rather than the 
first-principles part of the Hamiltonian often dominates the 
physics of the material. Attempts to calculate Hubbard U's 
or other parameters of the models are based on intuition and 
do not provide any solid foundation for these models since 
any connection between the calculations and the models are 
tenuous at best. As far as we are aware, the only reasonable 
first-principles method for calculating electronic correlation 
effects in metals beyond LDA is the GW approximation. Al- 
though this is a low-order approximation that definitely fails 
for very strong correlation effects, it is sufficient for our pur- 
poses as a way to estimate correlation deviations from LDA 
band-structure theory. The main purpose of our work is to 
show that it is possible and useful to define a new quantitiy, 
which we call correlation strength, for both experimental and 
theoretical work on new materials in order to place these sys- 
tems in their proper physics context. 

As a theoretical method for estimating correlation effects 
we have used the quasi-particle self-consistent GW approx- 
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imation (QSGW) EMI- 

The GW approximation can 
be viewed as the first term in the expansion of the non- 
local energy-dependent self-energy S(r, r', oj) in the screened 
Coulomb interaction W. From a more physical point of view 
it can also be interpreted as a dynamically screened Hartree- 
Fock approximation plus a Coulomb hole contribution Eli. 
Therefore, GW is a well defined perturbation theory. In 
its usual implemention, sometimes called the "one-shot" ap- 
proximation, it depends on the one-electron Green's functions 
which use LDA eigenvalues and eigenfunctions, and hence 
the results can depend on this choice. Unfortunately, as cor- 
relations become stronger serious practical and formal prob- 
lems can arise in this approximation [II 1Q. However, Kotani et 
al. II 1 211 have provided a way to surmount this difficulty, by us- 
ing a self-consistent one-electron Green's function that is de- 
rived from the self-energy (the quasi-particle eigenvalues and 
eigenfunctions) instead of LDA as the starting point. In the lit- 
erature, it has been demonstrated that the QSGW form of GW 



To set an appropriate correlation scale, we define our theo- 
retical C by: 



theory reliably describes a wide range of semiconductors 11311 . 
spd ifiol Q, fl5tl and rare-earth systems flill . It should be 
noted that the energy eigenvalues of the QSGW method are 
the same as the quasiparticle spectra of the GW method. This 
captures the many-body shifts in the quasiparticle energies. 
However, when presenting the quasiparticle DOS, this ignores 
the smearing by the imaginary part of the self-energy of the 
spectra due to quasiparticle lifetime effects, which should in- 
crease as quasiparticle energies become farther away from the 
Fermi energy. 

To define a theoretical correlation strength some electronic- 
structure quantity that scales with an intuitive notion of corre- 
lation strength is needed. In our application to Pu, we propose 
to consider the /-bandwidth, Wf, and use the relative band- 
width reduction in QSGW compared to LDA, 



w rei = W)(GW)/W>(LDA), 



(1) 



as the key quantity, where Wf(GW) and Wf (LDA) are the /- 
bandwidths as obtained from QSGW and LDA calculations, 
respectively. This is consistent with the correlation-induced 
QSGW /-bandwidth reduction in Pu that was demonstrated 
inRef.S 

Using a quasiparticle calculation is important since lifetime 
effects, which are absent in the LDA calculations would ob- 
scure the band narrowing in GW relative to LDA. We also 
need a measure that is robust at the high temperatures of the 
strongly correlated phases of Pu, where any low energy fea- 
tures in the electronic structure are likely to be thermally av- 
eraged away. II 1711 . In this regard, it should be noted that, al- 
though temperature certainly plays an important role in pre- 
dicting the correct equilibrium crystal structure, we believe 
that it is the resulting volume per atom of any Pu phase that 
determines the amount of correlation, since this is an elec- 
tronic property. In particular, we don't expect that the band- 
width predicted by our zero-temperature GW calculations will 
be sensitive to any temperature in the range set by the Pu solid 
phases. 



C = 1 - w re 



(2) 



which ranges from C = (no bandwidth reduction) in the 
LDA limit to C = 1 in the fully localized or atomic limit (the 
bandwidth becomes zero). 

As mentioned above, our test case for correlation is ele- 
mental Pu, an actinide metal, which exhibits large volume 
changes compared to predictions from band structure the- 
ory that are clearly due to correlation effects 1 1 8142211 . The 
large variation in volumes is controlled by the amount of very 
strong /-bonding, which is due to direct /-/ wave-function 
overlap. The /-bonding for many of the different phases is 
greatly reduced leading to anomalous volume expansions due 
to the narrowing of the /-bands that results from correlation 
effects [22[| . If no correlation were present, the /-bonds would 
have their full strength and a relatively small volume per atom 
for all phases would be accurately predicted by LDA band- 
structure methods. In the limit of extremely strong correlation 
the bands would have narrowed so much that the / electrons 
would be fully localized, and they would not contribute to the 
bonding. The volume per atom would then be much larger 
and close to that of Am, which has fully localized / electrons 
that do not extend outside the atomic core. 



Using the QSGW approximation we have calculated 12311 
the quasi-particle band structures of the fee, bec, simple cubic 
(sc), 7, and the pseudo-a phases of Pu as a function of volume. 
The pseudo-a is a two-atom per unit cell approximation 02411 
to the true a structure of Pu that preserves the approximate 
nearest-neighbor distances and other essential features needed 
for the electronic-structure. In this way we avoid performing 
an extremely large and expensive 16 atom-per-unit-cell cal- 
culation for the a-structure. We are unfortunately unable to 
present GW results for the /3-structure, which is even more 
complex than the a structure, since no pseudo-structure for 
this crystal structure is available and a QSGW calculation is 
presently not feasible for so many atoms per unit cell. 

To calculate the /-electron bandwidths from the /-electron 
projected density of states (DOS), D/(E), an algorithm is 
needed to determine the width of the main peak in this DOS. 
A simple first guess is to choose a rectangular DOS and to 
use a least-squares fit to the GW or LDA /-DOS to deter- 
mine the best height and width of the rectangle. A drawback 
of this method is that an artificial broadening of the effective 
/-bandwidth appears, which is due to a significant d-f hy- 
bridization at the bottom of the /-DOS that creates an extra 
peak at low energies. This masks the correlation induced band 
narrowing. Since this peak has relatively lower height than the 
main / peak, we may avoid this complication by generating an 
algorithm that emphasizes the "high-peak" part of the /-DOS. 
The algorithm we have used is therefore the second moment 
of the /-DOS 



W = 2((E 2 ) - (E) 2 ) 1 / 2 . 



(3) 



3 



P Y e 8 
I I I I 






5 




e 




sc 




ps-a 




Y 



14 16 18 20 22 24 26 28 



V (A 3 ) 



FIG. 1. (Color online) Plot of w rei = W f (GW)/W f (LDA) versus 
volume, V, per atom, for the 7, fee, bee, sc, and ps-a (pseudo-a, 
an approximate a-phase f24ll ) crystal phases of Pu. Note that the sc 
(simple cubic) is a hypothetical structure for Pu. The small, vertical 
bars at the top of the figure mark the experimentally observed atomic 
volumes [26]. 
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FIG. 2. (Color online) Plot of w rei = W/(GW)/VK/(LDA) versus 
W/(LDA) for the 7, fee, bee, sc, and ps-a. The dashed red line 
represents the fit of Eq. (O The small, vertical bars at the top of the 
figure mark the values of Wj (LDA) calculated at the experimental 
volumes of the five Pu phases [26]. 



The factor of two is needed because the bandwidth extends 
above and below the mean energy and is not just the average 
deviation from the mean energy. To emphasize the main part 
of the /-DOS peak, the square of the /-DOS is used as weight 
function I25I1 : 



i.e., 
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J dEf(E)D}(E)/ J dED%E). (4) 



In Fig.[T]we illustrate how w re ; varies with volume for the 
five different phases considered here lErjll . Large volume vari- 
ations ranging between about 14-28 A 3 per atom are consid- 
ered, with bandwidths that span almost an order of magnitude, 
from ?»0.5 eV to «2.5 eV. Although the LDA bandwidth de- 
creases with increased volume due to reduction in /-/ over- 
lap of the wavefunctions, the QSGW bandwidth decreases 
even faster illustrating increased correlation effects with lat- 
tice expansion. The bandwidth at a specific volume depends 
on crystal structure (due to differences in coordination and 
bond lengths), as does also the correlation strength. 

Although we expect electronic-structure calculations to 
strongly depend on the crystal structure and lattice constant, 
we surprisingly found that correlation effects were approxi- 
mately independent of these. Indeed, Fig. f2] shows that all 
of our different calculations for our measure of correlation 
strength, the reduced bandwidth, collapse to a single "uni- 
versal" curve when plotted as a function of the LDA band- 
width. In making this plot, it is likely that the effective 
screened Coulomb interaction between the 5/ electrons is ap- 
proximately constant and that the correlation effects are being 
tuned by the effective average kinetic energy of these elec- 
trons as reflected in their LDA bandwidth. In the range of Wt 
values considered here the curve is approximately quadratic, 



w, 



ei(x) = 0.15 + 0.43.x - 0.07a: 2 



(5) 



where x = W/(LDA) in eV. From Eq. (O we can use these 
results to determine a correlation strength C. It is remarkable 
that the many-body properties of a strongly correlated system 
can be tuned with what is normally considered to be a one- 
electron property. 

In Fig. [3] we show i27ll that our definition of theoretical 
correlation strength does indeed fulfill our expectations and 
can be used to bring order into the trends for various ex- 
perimental properties, including volume, sound velocity, and 
resistivity. These properties exhibit an approximately 25%, 
50%, and 35% change over the correlation range (about 0.2 to 
0.6) between the a and <5 phases of Pu and, with some scatter 
that might partially depend on sample quality, fall on smooth 
curves when plotted as a function of our theoretical correla- 
tion strength. It is remarkable that all of this data should col- 
lapse to a single curve for each property that is independent 
of any explicit consideration of temperature, crystal structure, 
or other variable. However, more generally, we would only 
expect this to be true for a property that was predominantly 
affected by correlation effects. 

In terms of theoretical trends, various theories have often 
attempted to estimate the amount of correlation in terms of 
the Z-factor, 



1 - (*, 



duj 



|*r,k) 



(6) 



where ^ n k are the (LDA) electronic eigenfunctions with en- 
ergies e„k, and E denotes the self-energy. We have found 
that the volume dependence of the Z-factors follows the trend 
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FIG. 3. (Color online) Trends in Pu properties as a function of cor- 
relation strength C, including (a) volume per atom [26], (b) sound 
velocity |2^1 . and (c) resistivity |2^1 . 



of the /-bandwidth reduction in Fig. Q] i.e., our measure 
of correlation strength, albeit with variations due to k- and 
hybridization-dependence. However, it should be noted that 
the relation between Z and bandwidth reduction is not the 
same in all materials, especially for weakly correlated broad- 
band systems, which seem very different from strongly corre- 
lated materials like Pu. 

The simplest Hubbard-like Hamiltonian to describe 
strongly correlated electron systems has a form 



H 



(7) 
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FIG. 4. (Color online) C from GW theory versus 1/W f(LDA). The 
data for Co, Rh and Ir are for the 3d, 4d, and 5d bandwidths, re- 
spectively. The small, vertical bars at the top of the figure mark the 
values of W/(LDA) _1 calculated at the experimental volumes of the 
five Pu phases [26]. 



with two parameters: the Hubbard parameter U which induces 
correlation, and an effective t, which can be related to the un- 
correlated bandwidth W . When W dominates, the system 
is in a weakly correlated limit and, when U dominates, the 
system is in a strongly correlated regime. Hence, one can 
study the solutions as a function of U/W to go from one 
limit to another. In more realistic electronic-structure calcu- 
lations, the same physics is intuitively expected to carry over. 
The Hubbard U can then be thought of as a screened onsite 
Coulomb interaction and the bandwidth as due to the normal 
band-structure hybridization. In our context, this suggests that 
the correlation strength C should also be a function of U /W. 
To test this, in Fig. |]we plot C versus 1/W/(LDA). If the 
effective U were approximately constant, we had hoped to ob- 
serve some approximate linear behavior at weak correlations, 
but any such behavior is unclear in Fig.|4] To show what might 
happen at weaker correlation strengths we have also included 
in Fig.|4]also the equilibrium-volume results for Co, Rh, and Ir 
for the d-electron projected DOS. Interestingly enough, the d- 
electron results seem to follow the same overall trend to large 
bandwidths (small correlation). Among the transition metals 
included in the plot, Co (3d) has the most narrow d-band, and 
the correlation value is close to the lowest values for Pu in the 
figure. 

In summary, we have introduced the idea of a "correla- 
tion strength" quantity C, which must be taken into account 
in order to explain the properties of strongly correlated elec- 
tronic materials. As an example, we have shown how to use 
the GW method to define a theoretical C for metallic Pu, 
and that various experimental physical properties, including 
anomalous volume expansion, sound velocity, and resistiv- 
ity, for the different phases of Pu follow well defined trends 
when plotted versus our theoretical correlation strength. We 
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have also demonstrated a universal scaling relationship for the 
correlation-reduced bandwidth as a function of the LDA band- 
width. 

This work was carried out under the auspices of the Na- 
tional Nuclear Security Administration of the U.S. Depart- 
ment of Energy at Los Alamos National Laboratory under 
Contract No. DE-AC52-06NA25396, the Los Alamos LDRD 
Program, and the Research Foundation of Aarhus University. 
The calculations were carried out at the Centre for Scien- 
tific Computing in Aarhus (CSC-AA), financed by the Danish 
Centre for Scientific Computing (DCSC) and the Faculty of 
Science and Technology, Aarhus University. 



* Electronic address: svane@phys.au.dk 
[1] L. Hedin, J. Phys.: Condens. Matter 11, R489 (1999). 
[2] F. Aryasetiawan and O. Gunnarsson, Rep. Prog. Phys. 61, 247 
(1998). 

[3] A. N. Chantis, R. C. Albers, A. Svane, and N. E. Christensen, 

Phil. Mag. 89, 1801 (2009). 
[4] A. Kutepov, K. Haule, S. Y. Savrasov, and G. Kotliar, Phys. 

Rev. B 85, 155129 (2012). 
[5] A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Rev. 

Mod. Phys. 68, 13 (1996). 
[6] G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, O. Par- 
collet, and C. A. Marianetti, Rev. Mod. Phys. 78, 865 (2006). 
[7] K. Held, O. K. Andersen, M. Feldbacher, A. Yamasaki, and 

Y.-F. Yang, J. Phys.: Condens. Matter 20, 064202 (2008). 
[8] J. Kunes, I. Leonov, M. Kollar, K. Byczuk, V. Anisimov, and 

D. Vollhardt, Eur. Phys. J. Special Topics 180, 5 (2010). 
[9] R. C. Albers, N. E. Christensen, and A. Svane, J. Phys.: Con- 
dens. Matter 21, 343201 (2009). 
[10] M. van Schilfgaarde, T. Kotani, and S. Faleev, Phys. Rev. Lett. 

96, 226402 (2006). 
[11] M. van Schilfgaarde, T. Kotani, and S. V. Faleev, Phys. Rev. B 

74, 245125 (2006). 
[12] T. Kotani, M. van Schilfgaarde, and S. V. Faleev, Phys. Rev. B 

76, 165106 (2007). 
[13] A. Svane, N. E. Christensen, M. Cardona, A. N. Chantis, M. van 

Schilfgaarde, and T. Kotani, Phys. Rev. B 81, 245120 (2010). 
[14] S. V. Faleev, M. van Schilfgaarde, and T. Kotani, Phys. Rev. 

Lett. 93, 126406 (2004). 
[15] A. N. Chantis, M. van Schilfgaarde, and T. Kotani, Phys. Rev. 

Lett. 96, 086405 (2006). 
[16] A. N. Chantis, M. van Schilfgaarde, and T. Kotani, Phys. Rev. 

B 76, 165126 (2007). 
[17] The most strongly correlated phase of Pu is 5 Pu, which 
has a temperature between about 600 and 700K. As noted in 
Ref. l28ll , "the derealization process of the 5f electrons," i.e., 
the correlation effects, continue to produce anomalies as high 



as 2000 K in temperature in many of these properties, well into 

the high- temperature liquid phase of Pu. Also, see Ref. l29Tl . 
[18] O. J. Wick, Plutonium Handbook: A Guide to the Technology 

(Gordon and Breach, New York, 1967). 
[19] S. S. Hecker, MRS Bull. 26, 672 (2001). 
[20] S. S. Hecker, Met. Mat. Trans. A 35A, 2207 (2004). 
[21] S. S. Hecker, D. R. Harbur, and T. G. Zocco, Prog. Mat. Sci. 

49, 429 (2004). 
[22] R. C. Albers, Nature 410, 759 (2001). 

[23] We have not included spin-orbit effects, which can be safely 
ignored for the purposes of this paper. The Pu /-DOS splits 
into a pair of clearly separated j=5/2 and 7/2 peaks. To include 
spin-orbit, we would need to calculate the bandwidth of each 
peak separately and use that corresponding to j=5/2. By ignor- 
ing spin-orbit coupling, we are saved from this additional trou- 
ble, which is not expected to change the effective /-bandwidths. 
Recent spin-orbit GW calculations have been calculated in Pu 
0]. However these have been done in the fully self-consistent 
GW method, which usually is a poor approximation in solids 
due to an incorrect treatment of plasmon effects. Since the DOS 
in this paper includes broadening effects due to the imaginary 
part of the self-energy in all of the different approximation that 
were used, it is also unclear how bandwidth narrowing would 
separately be affected by spin-orbit effects. 

[24] J. Bouchet, R. C. Albers, M. D. Jones, and G. Jomard, Phys. 
Rev. Lett. 92, 095503 (2004). 

[25] Choosing instead for example Df(E) as the weight function 
does not serve the purpose of emphasizing the central over the 
"hybridization wings" in the /-DOS. 

[26] For the atomic volumes we have ignored any thermal volume 
expansion. Each phase is represented by a volume correspond- 
ing to a fixed temperature within that phase. We have used the 
original data of Zachariasen and Ellinger [3G--33] correspond- 
ing to the volumes at the temperatures: 21, 190, 235, 320, 477, 
and 490 °C, for the a, /3, 7, S, 5' , and e phases, respectively. 

[27] For the volumes of the different phases of Pu, we have followed 
the same method used to generate Fig. [TJ We have also used 
the same volumes of the different phases for the sound velocity 
and resistivity needed to determine the correlation strength from 
the GW calculations plotted in Fig.Q] Note that, since we have 
not directly calculated the value of w re ; for the j3 phase, we 
instead used the availability of the bandwidth reduction of Eq. 
(O together with the calculated LDA bandwidth for the correct 
crystal structure of /3-Pu to determine w re i (/3) = 0.55. 

[28] M. Boivineau, J. Nuc. Mater. 297, 97 (2001). 

[29] M. Boivineau, J. Nuc. Mater. 392, 568 (2009). 

[30] W. H. Zachariasen and F. H. Ellinger, Acta Cryst. 16, 777 
(1963). 

[31] W. H. Zachariasen and F. H. Ellinger, Acta Cryst. 16, 395 
(1963). 

[32] W. H. Zachariasen and F. H. Ellinger, Acta Cryst. 8, 431 (1955). 
[33] F. H. Ellinger, J. of Metals 8, 1256 (1956). 



